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Abstract 

We introduce a novel method to find exact density operators for a spin-1/2 particle in time- 
dependent magnetic fields by using the one-mode bosonic representation of su{2) and the connection 
with a time-dependent oscillator. As illustrative examples, we apply the method to find the density 
operators for constant and/or oscillating magnetic fields, which turn out to be time-dependent in 
general. 
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I. INTRODUCTION 



Recently two 



lating fields 



evel atoms have been investigated under the influence of pulsed or oscil- 
However, a system of spin-1/2 particles under external magnetic fields 



has been a subject of interest since the advent of quantum mechanics p[. Also, two-level 
systems of electric dipoles under external electric fields have been studied for a long time 4 1 , 
and the statistics of a spin-1/2 particle has been found in terms of density matrices 5|. 
All these systems can be described by the Lie algebra su{2) whose generators satisfy the 
commutation relations for angular momentum (in units oih = 1) 

[Sa, Sb] = ieabcSc, (1-1) 

where eabc is the Levi-Civita tensor, taking +1(— 1) for even (odd) permutations of 1, 2, 3 
and vanishing otherwise. The Lie algebra su{2) can be used to find the evolution operator 
for the spin-1/2 particle system j^. 

On the other hand, a harmonic oscillator has the Lie algebra su(l, 1), which is isomorphic 
to so{2, 1) j?]. The generators of su{l, 1) produce the spectrum of a number of quantum 



systems. These generators satisfy the commutation relations 

[K,, K2] = -tK,, [K2, K,] = zK,, [K,, K,] = 1K2. (1.2) 



We may introduce the position representation of sm(1, 1) as 

*--K4^^)- f^-^' 

In the above definition, we scaled y = 2q as in Ref. i^J- The Lie algebra su{l, 1) has also 
been used to find the quantum states of an oscillator jsf. 

The compact Lie algebra su{2) is isomorphic to so(3), (5^(2) ^ so(3)), whereas the 
noncompact Lie algebra su{l, 1) is isomorphic to so(2, 1), (su(l, 1) ~ so(2, 1)) [7]. In fact, 
the compact group SO {3) generated by so(3) holds x\ + x\ + x\ invariant whereas the 
noncompact group 50(2, 1) generated by so(2, 1) holds x\ + x\ — x\ invariant. Though 
the spin generators of 5^(2) satisfy the commutation relations for an angular momentum 
operator, many years ago Schwinger introduced the two-mode bosonic representation for 
su{2) 0. It is, thus, tempting to exploit the connection between the spin particle and the 
oscillator. As su{2) is compact but sm(1, 1) is noncompact, the connection between these two 
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algebras requires a complexification of generators. This can be achieved by introducing a one- 
mode bosonic representation of su{2), Sx = —K3, Sy = —iKi, and Sz = iK2. The generators 
Sy and Sz are anti-Hermitian because Ki and K2 are Hermitian operators. Nevertheless, it 
is likely that the oscillator may be used to find the evolution of the spin particle just as the 
evolution operator of an oscillator can be found from a spin-1/2 system jg,]. 

The density operator of a physical system is an important tool for understanding the 
statistical properties of a system. The density operator, p(t), provides an ensemble average, 
Tr[p(t)A], for every physical observable A. In this paper, we introduce a novel method to 
solve the homogeneous Bloch equation from the connection with the time-dependent os- 
cillator that uses the one-mode bosonic representation of su{2) and find a class of density 
operators for a spin-1/2 particle in time-dependent magnetic fields. A time- dependent oscil- 
lator is known to possess invariant operators, or density operators, that satisfy the Liouville 
equation Q]. We use some known invariant operators for the oscillator to find the density 
operators for the spin particle, for instance, in constant and/or oscillating magnetic fields. 

The organization of this paper is as follows: In Sec. II, we briefiy discuss the density 
operator for a spin-1/2 particle in a time-dependent magnetic field. The density operator 
satisfying the Liouville equation is determined by a solution of the Bloch or Landau-Lifshitz 
equation. In Sec. Ill, we introduce a one-mode bosonic representation of su{2) and find 
the oscillator that corresponds to the spin-1/2 particle. Further, using an invariant operator 
for the oscillator, we find a density operator of the spin particle through a connection of 
Lie algebras sm(1, 1) and 5^(2). In Sec. IV, we compare our method with the evolution 
operator and the complex parameter method. In Sec. V, we apply the method to find 
density operators in constant and/or oscillating magnetic fields and compare them with the 
well-known results. 

II. SPIN-1/2 PARTICLE IN A MAGNETIC FIELD 

We consider a spin-1/2 particle under the infiuence of an external magnetic field. For the 
sake of simplicity, we set 7 = g^B = 1? absorbing it into the magnetic field, where g and 
are the Lande g-factor and the Bohr magneton, respectively, of the spin particle. Then, the 
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spin particle has the Hamiltonian 

H{t) = - Mt)Sa, (2.1) 

a=x,y,z 

where Sa is the ath spin component. The spin particle obeys the time-dependent Schrodinger 
equation for a spinor state: 

t^^\^,t) = nim,t). (2.2) 

It is well understood that any instantaneous eigenstate of a time-dependent Hamiltonian 
is not, in general, an exact quantum state. Nor does the evolution operator follow U{t) = 
exp[—i JTt{t)dt] because [H(t'),H{t)] 7^ for t' 7^ t. We may directly solve either the 
Heisenberg or Schrodinger equation. However, instead of solving the Heisenberg equation, 
in this paper, we use the invariant operator method introduced by Lewis and Riesenfeld Q| . 
In fact, an invariant or density operator satisfying the Liouville equation 

^^ + m,m] = o (2.3) 

leads to a general quantum state of the form 

=^c„e'/<«'*N^-^l".*)|n,t), (2.4) 

n 

where 

I{t)\n,t) = X^\n,t). (2.5) 

For the spin 1/2-particle with su{2), one may find a density operator of the form 11 1 

I{t) = - Yl Ma{t)S,, (2.6) 

a=x,y,z 

where the vector M with components Ma satisfies the Bloch or Landau-Lifshitz equation 

'^^+MxB = 0. (2.7) 



dt 

Note that is a constant of the motion. There have been some attempts to find the 
evolution operator or density operator for a spin particle in magnetic fields 
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i^Q- 
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III. CONNECTION WITH AN OSCILLATOR 



We now introduce a one-mode bosonic representation for su{2): 

S. = l{p' + q'), Sy='-{p'-q'), S.= '-{pq + qp), (3.1) 

where the standard commutation relation [q, p] = i holds. Note that our one-mode bosonic 
representation of su{2) consists of Sx = —K^, Sy = —iKi, and Sz = iK2. Therefore, the gen- 
erators Sy and Sz are anti-Hermitian because Ki and K2 are Hermitian operators whereas 
and are both Hermitian. The non-Hermitian nature of the one-mode representation 
follows from the fact that su{2) is compact whereas su{l, 1) is noncompact. Though the 
bosonic representation in Eq. p.ip indeed satisfies the spin angular momentum relation 
in Eq. (jl.lj) . the relation S"^ = J/4, however, does not hold for the bosonic representa- 
tion because the oscillator has an infinite number of states. The oscillator Hamiltonian 
corresponding to the spin Hamiltonian in Eq. 1)2.11) takes the form 

Hit)^B44 + Blit)l^'^«±^, (3.2) 



with 

1 
2 



B+{t) = 7^{Bx + iBy), B_ = 5;. (3.3) 



A passing remark is that the Hamiltonian in Eq. ()3.2j) is not a Hermitian operator because 
the Sa are constructed as non-Hermitian. This means that the spin directions of Sa do not 
pertain to a real physical space, but to directions in a complex space. Nevertheless, the most 
important point is that the algebraic relations hold whatever representation one may use. 
Any physical interpretation should given to real operators that are obtained from algebraic 
relations in the end. Though invariant under the parity operator, q —q and p —p, 
the Hamiltonian does not respect time-reversal symmetry, p —p and i —i, and thus 
cannot be an example for "PT-invariant quantum mechanics Q|. 

From now on, we shall work on the non-Hermitian oscillator in Eq. ()3.2|) . The spin 
particle or oscillator with B+ = is trivial, so this case will not be considered anymore. In 
;he case of 5+ 7^ 0, the oscillator in Eq. ()3.2j) may have an invariant operator of the form 
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a{t) 



u{t)p - iu{t) - i'^u{t) j q 



(3.4) 
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In fact, a{t) satisfies the Liouville equation 

.da{t) 



dt 



+ [a{t),H{t)] = 



(3.5) 



when M is a solution to the auxihary equation 

d f u 



dt\B, 



+ 



4 ^dt\ B+ 



u 

'bI 



(3.6) 



1/2 

or written as u = -B+ v satisfies the equation in canonical form 

2 



B2_ d UBJ2\ 3/i?+y l(B, 
4 ^dt\ B+ ] ^\B+) 2\B^ 



(3.7) 



A second independent solution to Eq. ()3.6p may lead to another invariant operator. However, 
the complex conjugate u* cannot, in general, be a solution to Eq. (j3.6|) due to the complex 
coefficient 5+. Thus, a) is not another invariant operator, in contrast with the Hermitian 
oscillator case. 

Our stratagem to find the density operators for the spin particle is first to find invariant 
operators which are quadratic in q and p and then to use the inverse relations in Eq. (9). 
For instance, may lead, using Eq. (9), to the density operator 



X{t) = -2 



Srr. 



1<U 



B 



+ 



B, 



. (3.8) 



Therefore, we obtain the density operator in Eq. ()2.6|) with the complex components 



f 9 fu — i^u\'^ 



B 



2iiu 



-Aiu 



+ / 



' ii — i^u 

B^ 



(3.9) 



A direct calculation shows that the components in Eq. ()3.9|) satisfy Eq. ()2.7|) for the density 
operator. Furthermore, the complex M. satisfies the constraint equation 



Mi + Ml + M\ 



0. 



(3.10) 



Note that Eq. ()3.10|) is a consistent condition because A^^ is a constant of the motion for 
Eq. ()2.7|) . so the eigenvalue of the invariant operator, here = M.^, should be constant. 
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The real part M(r) and the imaginary part M(j) of Ai individually satisfy Eq. ()2.7|) . a linear 
equation. We, thus, have two real density operators: 



(3.11) 



a=x,y,z a=x,y,z 

Note that the real and the imaginary parts have the same magnitude, but they are orthogonal 
to each other 



M 



(r) 



M 



0. 



(3.12) 



IV. COMPARISON WITH OTHER METHODS 



We shall compare the method in Sec. Ill with other methods: the evolution operator and 



the method recently introduced by Kobayashi 



16 



17| . The evolution operator is determined. 



for instance, by Eq. (6.23a) of Ref. |18j|. The Schrodinger equation leads the evolution 
operator U, 



U 



V 

iw* 



iw 
v*. 



v*v + w*w = 1, 



(4.1) 



to satisfy the equation 

dfx 
dtVBl 



B2 d fiBj2' 



[x 



V, w). 



(4.2) 



Note that Eq. ()4.2|) is the complex conjugate of Eq. ()3.6|) . thus implying x = u*. In this 
sense, Eq. ()3.6|) will determine not only the evolution operator but also the density matrix, 
as was shown in Sec. III. 

Kobaysahi recently introduced two methods for solving the Bloch equation. In his first 
method, the Bloch equation is transformed to the rotating reference frame This idea of 
a rotating frame is in essence similar to using the time- dependent creation and annihilation 



15 



3 



. These operators 



operators in the Fock space of time-dependent harmonic oscillators 
are chosen to satisfy the Liouville equation, Eq. ()3.5|) or Eq. ()2.3|) . so the eigenvalue of the 
number operator N{t) = d^{t)d{t) is a constant of the motion, and the exact quantum state 
in Eq. (j2.4j) is an instantaneous eigenstate of the number operator up to a time-dependent 
phase factor. 

On the other hand, in his second method Kobayashi introduced complex parameters for 



the magnetization [17[. As the magnitude |M| of the magnetization vector M is a constant 



of the motion, the magnetization vector can be normahzed: 

M 



m 



IMI 



ml + rriy + ml = 1. 



He then introduced two complex parameters: 



V 



rux + im. 



y 



1 - 

— iiTLy 



1 + m. 



rUx — ifriy 



nix + imiy 

1 + TTLz 



The Bloch equation, Eq. ()2.7p . in his notation reading as 



dm ^ 

— — h 7B X m = 0, 

at 



leads to the Riccati equation 



7, 



7, 



There is a similar Riccati equation for vj. In our notation, it is written as 
We linearize the Ricatti equation in Eq. ()4.7p by introducing a new variable 



75_ 



z .7r 

-D; 

2 



Finally, we can obtain 
d I 



dt V 7S_ 



+ 



d ( 1-13,12" 



0. 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



dt\ -fB_ 

Note that Eq. ()4.9p is the same as Eq. ()4.2|) for the evolution operator because we scaled 
7 as 7 = 1 in Eq. ()2.7|) and B ^ — B. This means that our method based on the connec- 
tion between a spin-1/2 particle in time-dependent magnetic fields and a time- dependent 
oscillator with time-dependent mass and/or frequency is equivalent to the complex param- 
eter method of Kobayashi. Hence, our method provides a group theoretical foundation for 
Kobayashi's method. 
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V. EXACT DENSITY OPERATORS IN CONSTANT AND/OR OSCILLATING 
MAGNETIC FIELDS 

We now apply our method to find the class of density operators in Eq. (j3.8|) for a spin 
particle in magnetic fields. To show the nontrivial nature of the density operator, we find 
the density operator in the trivial case of constant and/or oscillating magnetic fields and 
compare it with the known result. 



A. Constant Field 

First, we directly solve Eq. ()2.7p for the density operator. In the case of a constant field, 
M parallel to Bq satisfies Eq. (12. 7|) . The density operator with the components 

M, = MoBo,, My = MoBoy, = Mq^o,, (5.1) 

where Mq is an arbitrary constant, is proportional to the Hamiltonian itself. In the coordi- 
nate system with the z-direction along Bq, a second solution may be found: 

M^ = Mo±cosiBot + ip), My = -Mo±smiBot + ip), M, = Mq^, (5.2) 

where Mqj_, Mq^, and cp are constants that are determined by the initial conditions. The 
vector M rotates around the magnetic field due to the torque r = M x Bq. 

Second, using the connection with the oscillator, we may find the density operator given 
in Eq. ()3.9|) . We look for a solution to Eq. ()3.6|) of the form 

u^t) = fio+Mo/'e'(^°*+^)/' (5.3) 

for a positive constant Mq and a phase constant (f. The complex components of the density 
operator are given by 



[Bqx + iBQyY — {Bq — BqzY 



My = ^Mo[(Eo, + z5o,)2 + (5o-5o,)2 



^ 2 

= Mq{Bq., + iBQy){BQ - 5o,)e^(^«*+'^). (5.4) 



Finally, we obtain a real density operator 



3 

1{r){t)=-Y.^ir)aSa. (5.5) 
a=l 



where 



M, 



{r)x 



--Mo fio'x - BL - {Bo - Bo,f cos(5ot + ^) + MoBo^Boy sm{Bot + 



M, 



{r)y 



-Mn 



B^ - Bly + (5o - Bo,f sm{Bot + ^) - M^Bo^B^y cos{Bot + 



M, 



{r)z 



Mo{Bq - Boz) Bq^ cos(5ot + Boy sin(5ot + ip) 



(5.6) 



A few remarks are in order. Note that the components in Eq. (j5.6|) satisfy Eq. (j2.7|) or 
the Liouville equation. The two eigenvalues of the density operator in Eq. ()5.5|) are 



A± = T^M(,) = t\mMBo - Bo,] 



The imaginary components simply given by 



M 



{i) 



It 



-M 



(r) 



(5.7) 



(5.8) 



also satisfy Eq. ()2.7|) because the time derivative of a solution is another solution. Another 
solution to Eq. ()3.6|) . 

u_{t) = fio+Mo/'e-^(^«*+^)/2^ (5.9) 

leads to the components in Eq. ()5.(ij) for another density operator now with Bq replaced by 
-Bq and {BqI + ^) by -{B^t + 



B. Oscillating Field 

We turn to a spin particle in constant and oscillating magnetic fields 

B+ = Bo+e-'^\ B, = Boz, (5.10) 

where B^ = 2i?o+ cos(u;t) and By = — 2i?o+ sin(u;t). The oscillating magnetic field provides 
a resonance of the spin particle in the constant field Bo,- By directly solving Eq. ()2.7j) . we 
may find the density operator with the components 

l\/r 25o+Mo2 2Bq+Mqz . , X 

= , p cos(u;t), My = — — sm(u;t), = Mq^, (5.11) 

where Mq, is a constant. 

We now use the connection with the oscillator to find the density operator. Even the 
general case of B^ = -Bqx cos(co't) + Boysin^ut) and By = —BoxSin^ut) + BoyCos^ut) has a 
solution to Eq. ()H.fi|l of the form 



u+it) = 5o+Mo'/'e'[(^-")*+'^l/^ n = ^Bl + 2u;Boz + u;^. (5.12) 
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We then have the complex components 



Mx 



--Mo 
2 " 



-Mo 
2 



[Box + iBoy) e 



2 „—iujt 



{Box + iBoyfe-''^' + (^] - - 5, 



.2 Aujt 



M, = Mo{Box + iBoy){n - u - Bo,)e 



i{nt+ip) 



(5.13) 



From Eq. ()5.13j) . we finally obtain the density operator with real components 
1 



M, 



{r)x 



M, 



ir)y 



-Mo 
2 " 

--Mo 
2 ° 

--Mo 
2 " 

1 



|i?oa; — BQy — {Q, — uo — -Boz)^} cos{ujt) + 2BoxBoy sin(c<jt) cos(i7t + (p) 
ji^ox ~ BQy -j- {Q — uj — BqzY^ sin(c<jt) — 2BoxBoy cos{ujt) sin(fit + ip), 
[Bq^ - + {n-uj- Bozf} cos{ujt) + 2BoxBoy sin(cjt)] sm{nt + ip) 
+ ^Mo[{b^x -Boy + in-uj- Bo,f} sm{ut) - 2BoxBoy cos{ujt)] cos{Qt + <p), 



M(^r)z = Mo{fl — u) — Boz) Box cos(f2t + 9?) — Boy sin(fit + ip) 



(5.14) 



We, thus, find the general time-dependent density operator. Note that in the limit of = 
we recover the result in Eq. ()5.6|) for a constant field. Using an independent solution 



(5.15) 



we can find another density operator obtained by replacing Vt by —VL and [VLt + ip) by 

-int + ip). 

It would be interesting to compare the solutions Eqs. ()5.(ij) and ()5.14j) with those in 
Sec. Ill of Ref. ^]. However, our procedure to get solutions is simpler than that in Refs. 
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VI. CONCLUSION 



A spin-1/2 particle has a compact Lie algebra su(2), which can be connected through 
complexification with the oscillator algebra sm(1, 1). We introduced a non-Hermitian one- 
mode bosonic representation p.ljl for su{2). Using the bosonic representation, we found 
the oscillator in Eq. ()3.2|1 corresponding to the spin particle in an external magnetic field. 
The spin particle in the time-dependent magnetic field corresponds to a time- dependent 
oscillator. It is well known that the time-dependent oscillator has invariant operators which 
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provide exact quantum states up to time-dependent phase factors. A caveat is that the 
bosonic representation for su{2) is non-Hermitian and, as a consequence, the corresponding 
oscillator has, in general, a non-Hermitian Hamiltonian. However, this connection still 
provides us with a novel method for finding the nontrivial density operators for the spin-1/2 
particle, even for the well-known case of constant and/or oscillating magnetic fields. 

We have made use of the connection between the spin-1/2 particle and an oscillator to 
develop a novel method for finding the density operator in Eq. ()3.9p for the spin particle. The 
solution of the auxiliary equation, Eq. ()3.6|) . leading to the density operator is found to be 
the complex conjugate of the solution of Eq. ()4.2j) for the evolution operator, which in turn 
is the same as the solution of Eq. (j4.9|) for the complex parameters of the magnetization. 
As illustrative examples, we applied the method to the spin particle in constant and/or 
oscillating magnetic fields. The density operators in Eqs. ()5.6|) and ()5.14j) have a complicated 
time-dependence. In the case of the constant field, the density operator in Eq. ()5.1|) or 
()5.2|1 is widely used whereas the density operator in Eq. ()5.6|1 is the most general time- 
dependent one, which can also be obtained by using Kobayashi's method. Similarly, the 
density operator in Eq. ()5.14j) differs from that in Eq. (jS.llll and seems to be the most 
general one for the oscillating field. Thus, the connection with an oscillator provides an 
effective method for finding the general density operator for a spin-1/2 particle. 

The oscillator with a constant mass and frequency is known to possess not only time- 
independent invariant operators but also time-dependent ones, which lead to squeezed 
states of time- independent states [3]. The physical meanings and applications of the time- 
dependent density operator for constant field and/or oscillating magnetic fields will be ad- 
dressed in a future work. The method of this paper can be easily generalized to spin chain 
systems. Each spin particle corresponds to an oscillator; thus, a spin system is equivalent 
to a system of coupled oscillators. The invariant operator for the oscillator chain with time- 
dependent masses, frequencies, and couplings may be found in a similar manner 2(J]. It 
would be interesting to study the spin chain system in connection with the oscillator chain 
system. 
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